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We give a short personally-biased review on the recent progress in our understanding 
of gravitational radiation reaction acting on a point particle orbiting a black hole. The 
main motivation of this study is to obtain sufficiently precise gravitational waveforms from 
inspiraling binary compact stars with a large mass ratio. For this purpose, various new 
concepts and techniques have been developed to compute the orbital evolution taking into 
account the gravitational self-force. Combining these ideas with a few supplementary new 
ideas, we try to outline a path to our goal here. 

§1. Introduction 

Ground-based interferometric gravitational wave detectors have started opera- 
tion and are also in the phase of rapid improvement. R&D studies of a space- 
based gravitational wave observatory project, the Laser Interferometer Space An- 
tenna (LISA),^) which observes in the mHz-band, are in rapid progress. There is 
also proposal for a DECi hertz Interferometer Gravitational wave Observatory (DE- 
CIGO/BBO),^)'^) which will be a laser interferometer gravitational wave antenna in 
space sensitive at / ~ 0.1 Hz. 

Among the promising targets for space interferometers, the most precise theo- 
retical prediction and observational measurement of gravitational waveform are ex- 
pected from the inspiral stage of binary systems comprising of a supermassive black 
hole {M ~ 10^~^Mq) and a compact object of solar mass (/x ~ 1 — IOMq).^^ Here 
the meaning of "precise" is two-fold. One is that the number of cycles is large and 
the other is that higher order post-Newtonian corrections are necessary. These grav- 
itational wave sources can provide the first high-precision test of general relativity 
in very strong gravitating regimes.^-* 

On the theoretical side, to obtain sufficiently precise templates, the standard 
post-Newtonian approximation^'') seems to be insufficient since the accessible highest 
order of expansion is practically limited. However, there is a natural expansion 
parameter in this system, that is the mass ratio ^/M. Therefore, in place of the 
standard post-Newtonian expansion, there is a possibility to develop perturbation 
theory based on the well-understood black hole linear perturbation theory^^^^^^^ for 
binaries with extreme mass ratios. To develop a method to compute waveforms on 
this line is our purpose of studying gravitational radiation reaction. 

This paper is organized as follows. In Sec. |^ we briefly discuss the very basics 
of black hole perturbation approach, followed by a discussion on balance argument 
about radiation reaction, given in Sec. 3. In Sees. 4 and 5, we consider radiation 
reaction in adiabatic approximation. The basic idea is explained in Sec. 4, while the 
recent developments, which have led to a very concise formula for the change rate 
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of Carter constant, are presented in Sec. 5. A method for long time integration is 
also discussed. In Sees. 6 and 7, we discuss the instantaneous self-force. The main 
stream of the procedure about how to evaluate the self-force, although biased by my 
personal view, is explained in Sec. 6, supplemented with comments on some topics 
given in Sec. 7. The discussions presented in Sec. 5 might be beyond the level of an 
introductory review paper. Section 8 is devoted to Summary. 
In this paper we use the units in which SttGn = 1- 

§2. Before radiation reaction 

We model a binary system by a point particle of mass fj, orbiting a black hole of 
mass M assuming fi <C M. In the lowest order in the mass ratio, 0((/i/M)*^), the 
particle moves along a geodesic on the background geometry. Using the black hole 
perturbation, we can compute the energy and angular momentum flux carried by 
gravitation waves emitted by the particle to the infinity or into the horizon. Already 
in this lowest order approximation, this black hole perturbation approach has proven 
to be very powerful for evaluating general relativistic gravitational waveforms. 

The perturbed metric is expanded as 

9^J.l^ — dfiu ~^ ^iiu ~^ ^\ii> + ■ ■ ■ • 
Then the linearized perturbed Einstein equations 

become coupled equations for metric perturbations h^-''^ in general, and are difficult 
to solve. However, when the background spacetime is given by Schwarzschild or 
Kerr black hole, perturbation equations can be written in terms of a single master 
equation due to symmetry. Here we schematically write the master equation as 

where L is a second order differential operator and C,^^^ is a master variable in the 
linear order. The source T^^^ is a function obtained from the energy momentum 
tensor. 

This reduction is called Rcggc-Whcclcr-Zerilli formalism in the Schwarzschild 
case,^^-*' and Teukolsky-Sasaki-Nakamura formalism in the Kerr case.^^-*'^^-* In 
the latter case, starting with metric perturbations is not very successful. Instead, 
we use Newman-Penrose quantities defined by, say, 

-2V' ~ -C^^pan^'m''nPfff, (2-1) 

as a master variable. This variable is a contraction of Weyl tensor Cfj,upa with two 
null tetrad bases, n'' (in-going) and ?7t^ (angular). Here note that is complex 
valued and " ~ " represents complex conjugation. We can similarly define 2V' by 
considering the contraction with i^^ (out-going) and m^. As a result, we have two 
equations 

,L,i; = ^sT, is = ±2) (2-2) 
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but they are not mutually independent. The source gT is obtained by applying a 
corresponding second order differential operator ^r^i/ to the energy momentum tensor 

Thus obtained master equation is separable. We introduce a harmonic function 
which describes angular dependence Sa{0, <^), where A = {i, m, u} is a set of two an- 
gular eigenvalues and one frequency eigenvalue. Homogeneous solution can be found 
in the form of a linear combination of mode functions s^A = sRAir) s5/i(^, if) e"**^*. 
Substituting this form, the problem to find a homogeneous solution reduces to solving 
a second order ordinary differential equation 

[d^. + ■ ■ ■] sRA{r) = 0. (2-3) 

By using homogeneous solutions which satisfy appropriate boundary conditions, 
we can construct Green function. Then, the solution with source is given by 

si^^Yl l^^'^r(^) / V^d''x's^T\x'mx')e{r -/) + •••, (2-4) 

where up and out, respectively, mean that there is no in-coming wave from the past 

infinity and that there is no wave absorbed into the future horizon. Here, "• • • " stands 
for the part which contains 9{r'—r). Wc have introduced spin-flipped mode functions, 
sf2A{x') = _,i?^(r)3.S^(^,^)e--*. Notice that = sR'^^iriTSAO^^)^. 

(Since the differential operator L for the Teukolsky equation is not real, a simple 
complex conjugation does not give a solution of the same equation but it becomes 
a solution of the equation with its spin flipped from s to —s. Therefore we cannot 
choose sR^^ = sR"^^ as usual, but we can choose s-R*" = -sR"^^- ) sWa is the 
Wronskian between these two radial functions, i.e., 

,Wa = W{sRj,sR'X{r)) « sR7ir)^^^^ - sR^irf-^^^- (2-5) 

There is a systematic method to solve homogeneous Teukolsky equation as a series 
expansion.^^)'^®)'^^) Using this method, one can rather easily obtain the solution, 
and we can write down the expression analytically in explicit form once we invoke 
the post-Newtonian expansion. 

To estimate the fluxes at r ^ oo, one can use a simple relation, _2'0 ~ \{^+ ~ 
ihx)- Using this relation, one can estimate, say, the energy loss rate due to gravita- 
tional wave emission per unit steradian as 

In this paper, we use £ {C) as the specific energy (angular momentum) per unit mass 
of the particle. The contribution from the waves absorbed by the black hole can be 
evaluated in a similar manner. 



§3. Balance argument and its limitation 



When radiation reaction is not very significant, the trajectory of a particle is 
almost a geodesic. In such a case one can estimate the evolution of orbital frequency 
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by assuming that the orbit evolves losing its energy and angular momentum as much 
as those emitted as gravitational waves. ^^-^ 

In the case of circular orbits, the evolution of the orbital frequency almost deter- 
mines the gravitational waveform (except for the change of amplitude) . Equating the 
energy lost through the gravitational wave emission with the minus of the binding 
energy, the change rate of the orbital frequency / is evaluated as 

df _ dSpw ( d£orbit \ ^ /„,N 
dt~ dt \ df ) ■ ' 

Then the leading order of dScw/dt starts with 0(/u). (Recall that £ is the specific 
energy. In counting of the order of magnitude, we are assuming that M is 0{\).) 
This is because the energy was exactly conserved, if the radiation reaction were 
turned off. On the other hand, the leading order of the relation between the binding 
energy and the orbital frequency, d£orbit/df, is determined by the geodesic motion 
on a given background black hole spacetime. The effect of the self-force, which is 
the force acting on a particle due to the metric perturbation caused by the particle 
itself, is higher order correction of 0(/x) in d£orbit/df. Roughly speaking, to obtain 
the leading order correction to the waveform, we therefore have only to know the 
effect of the self-force on d£Gw/dt. (We will return to this point in Sec. 

The extension to more general orbits is straight forward in the case of Schwarzschild 
background. In this case one can assume that the orbit is on the equatorial plane 
without loss of generality. Then, the geodesies are specified by the energy £ and 
the z-component of the angular momentum C. Both "constants of motion" have 
the associated timelike and rotational Killing vectors, r]^^^^ = [dt)^ and rj^^^ = {d^,)^. 
Therefore we can define conserved current from the effective energy momentum ten- 
sor t^^, which satisfies the conservation law t'^'^-i, = with respect to the background 
covariant derivative, as 

3f = [^-%vrr(^- (3-2) 

Then satisfies j^^^ = 0, and hence one can define the conserved (specific) 
energy by the integral over spatial three surface E as 

£ = - 

Let's choose the boundary of i7 to be a sphere S at a fixed radius, which is supposed 
to be sufficiently large. In this case the change rate of £ is given by 

^-j/S,P, (3.4) 

where dS is an infinitesimal element of area on S. In the above discussion, we 
assumed the existence of the effective energy momentum tensor t^j^ that satisfies 
the conservation law. Now we directly derive it from the Einstein equation. The 
Einstein tensor can be expanded as 

G^,[g + h] = Gfl + gW [h] + Ggj,[h, h] + • • • . (3-5) 




Gravitational radiation reaction 



5 



Since the background metric g satisfies the vacuum Einstein equations in the present 
case, g[S/ = 0. Using this relation, the contracted Bianchi identity at the hnear order 
in h becomes 

Gil[hr = (3-6) 

Here a semicolon denotes a covariant differentiation with respect to the background 
metric. Substituting "h = h^^^ + h^'^^ + • • • " into the Einstein equations and keeping 
the terms up to second order, we obtain 

M^)] + Gil[h('^] + Ggl [h(^),h(^)] = T,.. (3-7) 

From the background covariant derivative of this equation, we have 

Ggl[h«,h«]^'^ = r^/. (3-8) 

Hence we find that 

t^.^T^.-G^^l[hW,hW] (3-9) 

satisfies the conservation law with respect to the background covariant derivative. 
One may identify the second term as the effective energy momentum tensor of grav- 
itational waves t^iil^\ Accordingly, one can divide the energy £ into two parts: 
£ = £orbit + ^GW ■ ^orbit is the Contribution from T^^^ while £gw is that from t^^^^ . 

Now we should note that £orbit is 0(1) while £gw is 0{fi). Keeping this fact 
in mind, we look at Eqs. ()3-3() and 1)3 ■4() again. Then £ is 0(1) and the leading 
term comes from £orbit- On the other hand, d£/dt does not have contribution from 
T^y, and hence it is 0{fi). However, we do not conclude that the radiation reaction 
is unimportant. After integration over a long period of 0{fi~^) the change in £ 
becomes 0(1). This change in £ must be attributed to the change of T^^, because 
£gw should stay 0{fj,). This consideration establishes the balance argument to the 
motion of a particle: 



d£ 



orbit 



dt 



dSjP. (3-10) 



s 



In the case of Kerr background, we do not have spherical symmetry. Hence, one 
cannot say that the orbits are on the equatorial plane in general. To specify geodesies 
off the equatorial plane, we need to consider another "constant of motion" . It is well 
known that the geodesies on the Kerr background possess the third "constant of 
motion" called Carter constant. However, Carter constant is not associated with 
any Killing vector field. Instead, it is related to a rank two Killing tensor. 

Let us examine the "constants of motion" in Kerr in more detail. The back- 
ground Kerr spacetime in the Boyer-Lindquist coordinates is given by 

ds=-{l — \ dt^ — dtdLp 

, 2 , ^^a2 , 2 , 2 , 2Ma'^rWe \ . 2 2 .onx 
+—dr + Ld8 + lr+a-\ — I sm Odip , (3T1) 
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where Kerr parameter a is the the angular momentum of the black hole divided 
by its mass, i7 := + cos^ 9, and Z\ := — 2Mr + a^. Killing tensor in Kerr 
spacetime is given by 

K^^ := 217/(^71^) + r'^g^u, (3-12) 

where the parentheses denote symmetrization on the indices enclosed, and := 
(r^ + a^, Z\, 0, a) /A and := (r^ + a^, — Z\, 0, a) /2X! are out-going and in-going 
radial null vectors, respectively. Killing tensor satisfies the equation 

i^(M.;p) = 0. (3-13) 

Using this Killing tensor, the Carter constant is defined as 

Q = K^pu^u^, (3-14) 

where := dz^/dr is the four velocity of an orbiting particle. We often use another 
notation for the Carter constant, C := Q — {a£ — defined in such a way that it 
vanishes for orbits on the equatorial plane. By using the symmetry of the Killing 
tensor, it is easy to check that thus defined Carter constant does not vary along 
geodesic. For Carter constant, however, we cannot define a quantity corresponding 
to £gw- Therefore there is no counter part of Eq. (|3-10() for dQ/dt. 

§4. Adiabatic approximation for dQ/dt 

As we have seen in the previous section, one cannot use the balance argument 
to evaluate the change rate of Q. Then, to evaluate dQ/dt, we need to compute 
the self-force acting on the particle directly. ^'^^ Though the prescription to calculate 
the self-force is formally established, ■^^•'"^'^^ performing explicit calculation is not so 
straight forward. However, it was found to be easier to compute the averaged value 
of dQ/dt. As an approximation, we may use the averaged values of the change rates 
for the "constants of motion" instead of those evaluated by using the instantaneous 
self-force. We call it adiabatic approximation. This approximation will be as good as 
an estimate using the balance argument, and moreover is also applicable to dQ/dt. 

4.1. use of radiative field 

Gal'tsov'^^-* advocated to use the radiative part of the metric perturbation, which 
was introduced earlier by Dirac,^^) to calculate d£/dt and dC/dt. The radiative 
field is defined by half retarded field minus half advanced one. The divergent part 
contained in the retarded field is common to that contained in the advanced field. 
Therefore the combination of the radiative field is free from divergence. Hence as 
far as we discuss the self-force composed of the radiative field, we do not have to 
worry about how to regularize the force. The question however remains whether the 
radiative self-force correctly reproduces the result which is obtained by using the 
retarded field with the aid of appropriate regularization procedure. 

Gal'tsov has shown that the radiative field correctly reproduces the results ob- 
tained by using the balance argument for d£ / dt and dC / dt when they are averaged 
over an infinitely long time interval assuming a geodesic motion as a source of metric 
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perturbation. When /i is thought to be an infinitesimal expansion parameter, the 
trajectory of the particle follows exactly a geodesic at the lowest order. In this sense, 
in the lowest order, this averaging over a long period of time assuming a geodesic 
motion is justified when we evaluate the self-force to the lowest order in 

However, there had been no justification for applying the same scheme to dQ/dt 
until very recently. The breakthrough was brought by Mino, who gave a justification 
for applying the same scheme to dQ/dt. "^^^ (See also Ref. 27)). Namely, he has proven 
that 



d ^\ 1 ,. I r , dQ 



— Q ) = - lim — / dT-^F"" 

where F" [h^^'^'^^] is the self- force evaluated by using the radiative field h^'""'^) = 

The key observation in his proof is the invariance of the geodesies under the 
transformation: 

a-^-a, {t,r,9,ip) ^ {-t,r,9,-ip). 

All geodesies, unless the values of (<?,£, Q) are fine tuned, transform into themselves 
under the above transformation if we choose the origin of t and (p appropriately. 
Once we notice this symmetry, it is easy to show that 

(ret)\ I / rin\ {o-dv)\ 

) , (4-2) 





where the superscripts (ret) and {adv) indicate that the retarded field h^'"'^*) and 
the advanced field h^"'^") are used instead of h^'""'^) in Eq. ()4-l() . The relation H4-2() 
justifies the use of the formula (|4-1|) . 

The radiative Green function has a simple structure which does not contain any 
step function 9{r — r'): 

G^'''"^^{x x') = V ^- 

^ ' ' ^W{sRY,sR']i)W{sR''r'',sRT) 

X (^{sR'x,sRT)s^T''{^)s^T''{A + w{sRT'',sR7)s^T{x)shT^ , 

(4-3) 

To show this, let us start with the following defining expression of the radiative 
Green function for r > r'; 



w(.flT,.-R*i') M'(,if'°«-»,,fl5»') 



Since our goal is to obtain an expression in terms of down-field and out-held, we 
want to eliminate gR^i^) ^'^d sR°a^ {r'){= _s^^(r')) in Eq. (|4-4|) . Hence we expand 
and as 



out I n jydown 



= a + (3 sR 
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= 7 sR""" + S sR"""- (4-5) 

Taking the Wronskian of both sides of Eqs. (|4-5|) with appropriate radial functions, 
one can easily obtain 

Substituting these relations, the expression ()4-4() reduces to ()4-3|) . We can do an 
analogous reduction for r < r', and the result turns out to be the same as that 
for r > r' . Namely, the step functions which was present in the retarded and the 
advanced Green functions do not appear in the radiative Green function. This means 
that the radiative field is a source-free homogeneous solution. 

4.2. metric reconstruction^^ ^ 

In order to use the formula (|4-1|1 . we need to know how to reproduce the metric 
perturbations from the master variable ^V- the present case, what we have to deal 
with is a source-free homogeneous solution. This fact simplifies the reconstruction 
significantly *\ The method was originally given by Chrzanowski.^^^''^^^ Here, we 
present the basic idea of the derivation of reconstruction formula, neglecting details. 

We formally write the metric perturbation induced by a source T"^ as 

V(^) = / V^gd'x'G%%{x,x')T»P{x'). (4-6) 

We assume that T"l^{x') is localized within r' < rg. And we assume a factorized 
form of the tensor Green function as 

= E -j^^ {KlAx)nTMx')9{r -/) + •••), (4-7) 

where n^^, is the mode function for the metric perturbation, whose explicit form is 
unknown at this point. Since we do not know how to normalize the mode function 
for metric perturbation, we have introduced a constant ^AA to take care of this 
normalization. is the Wronskian for the corresponding mode function for the 

master variable defined by (|2-5|) . 

The master variable can be computed from the metric perturbation by applying 
a second order differential operator as s^A = sD^" Uahu- Hence, we have 

= = E Jl^/^Ai^) I V^d'x'nTUx')T^^{x'), (4-8) 

for r > ro, outside the source distribution. On the other hand, we can evaluate the 
same quantity by using the Green function for the master variable, which leads to 

si> = sD^'^h^.ix) = ^ -l^,f}J{x) f ^d\' shT\x'),T{x'). (4-9) 

^ s ^ A J 

*■* When we reconstruct the metric perturbation from a solution of the master variable within 
the source distribution, more delicate treatment is necessary.'"'-' Below we use the assumption of 
the factorized form of the tensor Green function. This assumption is not necessary if we follow the 
derivation given by Wald,^^-* and actually this assumption itself is not correct."^"' 
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Comparing these two expressions, we find 

J ^gd^x' 7I^^t(x')T'^'^(x') = sM j V^d'x' ,^77*(x').r^.T'^"(x'). (4-10) 

The same relation holds for the modes with other boundary conditions. Since this 
relation holds for arbitrary T^i, as long as it satisfies the conservation law, T'^'^-,j = 0, 
one can establish the relation 

HA^^ix) = ,Mst;, s^Aix'), (4-11) 

where sT^^u is the second order differential operator that is obtained by integration 
by parts from sf^i/. This equation still has an ambiguity of adding pure gauge term 
since its contraction with T'^^ vanishes after integration over four volume. 
The above is a very crude explanation about the reason why we can reconstruct 
the metric perturbation form the master variable just by operating a second order 
differential operator. More rigorous derivation was given by Wald.29) With the aid 
of Starobinsky-Teukolsky identity*^ we can explicitly show that 

s^A = sD^'^nA^.u = sD'^'sMsT;, sf^A, (4-12) 

simultaneously fixing the normalization constant ^M. 

When the source is composed of a point particle, i.e., when the energy- momentum 
tensor takes the form T^^'^ = J dT{—g)^^^'^ u^u'^6^{x — z{t)), the radiative metric 
perturbation is given by 

^(w{,B^}^,,Rr)sUX^J^{x) I ^0T(z(r)) 

where 

^7 = Eu'^u'^nf;",- = ,meu^u\t;, s^a% (4-i4) 

and (j)'^^ is defined in a similar manner. For future convenience, instead of the four 
velocity u^, we used u^, an extension of to a vector field, whose definition is 
given below Eq. (|5-4j) . The factor S{= + cos^ ^) is also introduced for future 
convenience. 

We can easily confirm that gD^'^h^^u'^^ with the substitution of the above ex- 
pression reproduces the same gip^^'^'^^ that is obtained by using Eq. (|4-4|) neglecting 
the last complex conjugate term. This term is necessary to make the expression 
real. To show that adding this term does not disturb the property of reproducing 
a more detailed discussion is necessary. ^ However, we will not go into such 
a technical detail here. 

*-* "sD'"' sT^y" reduces to a forth order differential operator which transforms the radial function 
-sRa to the spin-flipped one sRa- This identity is called Starobinsky-Teukolsky identity. 
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§5. Simplified dQ/dt formula 

Now we know how to compute dQ/dt in principle. However, actual implementa- 
tion of calculation is not so straight forward. In this section we introduce a simpler 
expression for the adiabatic evolution of Carter constant. ^■^^ 

5.1. property of geodesies in Kerr 

We first discuss in a little more detail about geodesies in Kerr spacetime: z"(r) = 
(t2(r), r^(r), 0^(r), Here r is the proper time along the orbit. Using a new 

parameter A defined by d\ := dr/Z", which was recently reintroduced by Mino in 
the context of radiation reaction, the geodesic equations are written as 

drz \ ^ „ / N f d cos ^"^ 



dxj =®("°'^^)' (^-^^ 

^ = -a(a^: sin^ e.-C) + ^i±^P(r.), '^-^ = -a£ + + ^P(r.),(5.2) 

dX A dX sm 9^ A 

where P(r) = £{r'^ + a^) - aC, fi(r) = [P(r)]2 - A[r^ + Q] and 0(cos6') = C - {C + 
a^(l - £'^) + C'^) cos^ 6 + 0^(1 - S"^) cos^ 9. It should be noted that the equation for 
the r-component and the one for the ^-component are decoupled when we use A. 
The solutions of the first two equations (|5-1() are periodic. We denote the periods by 
21^1 fir and 2t: /Qq^ respectively. The other two equations (|5-2j) are integrated as 

MA) = t«(A)+t(^)(A) + (^)A, 

c^,(A) = c^M(A) + ^W(A) + ^^^A, (5-3) 

where (• • • ) means time average along the geodesic. t^^'\X) := / dX{{r1 + a?)P{Tz) / A 
-{{rl + a?)P{T,)/A)} and t^^\X) := - J dX{a{a£ sin^ 9, - C)-{a{aS sin^ 9^ - £))}, 
are periodic functions with periods 27r/i7j. and 27r/i7e, respectively. Functions ip^^'^ 
and (/j^^^ are also defined in a similar way. 

5.2. simplified formula 

We start to simplify the expression of the formula for dQ/dt. The self- force 
which is defined by u'^u^-u = f^, is derived from the geodesic equation on a 
perturbed spacetime. Then f^ is basically given by —SF^pau'^u'^, where SF^pa is the 
contribution to the Christoffel symbol from the metric perturbation h. Taking into 
account the redefinition of the proper time so that u^Up^ = — 1 is maintained, we 
obtain 

/'^[h] := -1(5^'^ + u^'u^){Kp,„ + K^.p - hp„,,)uPu\ 
Using this force, the evolution of Carter constant is given by 
dQ 
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lim 2 



(5-4) 



where ^{x) = Su^vyh^i, /2 and {ut,Ur,U0,Uip) := {—£, ±-\/R{r)/A, ±y^0{cos6) 
/sm9,C). This vector field tt^ is an extension of the four velocity of a particle in 
the sense that it satisfies u^{z{X)) = u^{X), but it differs from the parallel transport 
of the four velocity. Using the fact that Ur and ug, respectively, depend only on r 
and 6, we can easily verify the relation, Ua-js = U/3-a- 

If we take the long-time average of Eq. (|5-4|l . the second term in the last line 
vanishes because it is a total derivative. Furthermore we can show that the long- 
time average of the third term also becomes higher order in fi by using the relations 
Ua;i3 = Uf3-a and -/^(^jyjp) • Finally, we obtain 

We can derive an analogous expression more easily for the energy loss rate as^^^ 



lim 

dX / T^oo 2T 



T^oo 2T 



dXU{-r]%.)fa[hf,u 

-T 

1 

lim ^z:^ I dX 



7?f,))a,^(x) (5-6) 



where t?^^ is the timelike Killing vector. 

In any cases, all we need to know is ^{x) for the radiative field. Using the 
formulas (|¥T^ and , we find 

'^^'^'^^-) = ^j |^E^^>iX(-) / dX^^SJA^) + ■ ■ ■ , (5-7) 

where the normalization constant A^*" is defined by A^*" := —2TTujiW {sR'^'^'^ , sR^) 
sM-^W{sRj,sR'A)~^W{sRf'"'^,sRT)~^- Hereafter we neglect the contribution 
from waves absorbed by black hole, since the extension is trivial. 

We will not repeat here the technical issues for further reduction of the formula 
discussed in Ref. 32). Instead, we just briefly mention essential points. First point 
is that the frequency to is discretized as 

/ dA>f)(z(AO) = 2^'^ ('^ - ^i' (5-8) 

n,.,ng 

with 

w;;-"" := {dt,/dX)~^ (m {d^z/dX) + UrQr + ne^^e) ■ (5-9) 
Here A represents a set of eigenvalues {i,m,nr,ng}. For (dS/dt), we obtain 

f) = -Ei^iP. ("«) 
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and similarly, 



dC 



_ \ — _ I 7 



m 



(5-11) 



In the above, we have set 2A^*" = 1 by rescaling the amplitude of radial functions 
sRa appropriately. Although it is not manifest from our definition of A^*", we can 
show that it is real. For {dQ/dt), writing down the expression in Eq. H5-5|) explicitly, 
we arrive at 



dX 



x=z(X) 



dX 



-^{{r^ + a^)dt + ad^) - ^5,) (x) 



x=z{X) 



(5-12) 



Using Eqs. (|5-7| ) and (|5-8|) . the above expression reduces to 



3f? 



dX 



Pjr) 
A 



dX 



(5-13) 



x=z(X) 



The second point is to notice that the integrand is now a double periodic function 
with periods 27r/i7r and 27r/f?g. Suppose / is a function of g^^^ and g^^\ and 5^'"'* (A) 
and g^^\X) are periodic with periods 211 jQ^ and 2tt/Qq^ respectively. Then, in 
general. 



1 



hm — / dXf{g(^-\^),9^'\^)) 

dXr / dXef{g^'-\Xr),g^'\Xe)), 
10 Jo 



(2^)2 



(5-14) 



holds. Using this formula, the integral ()5-13|) can be integrated by parts. Then 
finally we arrive at the formula 



dQ 
dt 



(r2 + a^)P 
A 



d£\ _ / aP\ / dC\ ^sr^ n^Qr 



A 



dt 



This expression is as easy to evaluate as {d£/dt) and {dC/dt). To evaluate the last 
term, we have only to replace m with Urf^r in the expression for (dC/dt), ()5-ll() . 

Although the final expressions for (dC/dt) and {dQ/dt) are quite similar, there 
is a big difference between them. The squared amplitude of each partial wave, 
is a measurable quantity in the asymptotic regions, i.e., near the future null infinity 
or the future event horizon. The eigenvalues uj and m also can be read from the 
waveform in the asymptotic region. Hence, the expressions for {d£ /dt) and {dC/dt) 
are solely written in terms of the asymptotic waveform. This is consistent with the 
fact that the balance argument applies for these quantities. On the other hand, the 
eigenvalue iir, which appears in the expression for {dQ/dt), is not a quantity which 
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can be read from the asymptotic waveform without knowing the particle orbit. We 
can read the frequencies from the asymptotic waveform, but the frequency itself does 
not tell the numbers and ng, without an input of additional information about the 
orbit. This is consistent with our understanding that the balance argument cannot 
be used for evaluating (dQ/dt). 

5.3. long period orbital evolution 

One may ask how we can use the knowledge about the adiabatic evolution of Q 
to evaluate the long time orbital evolution. Just to be consistent with a given set 
of constants of motion we can choose r, 9 and ip arbitrarily at each time t as far 
as r and 6 are within the allowed range. Hence, the evolution of the constants of 
motion is not complete at all as a description of the orbital evolution. Nevertheless, 
combined with the normalization of four velocity, 

nrr dz^^ dz^ 

three constants of motion are sufficient to specify the four velocity as a function of 
r and 9. Even if we take into account the instantaneous self- force, the equations of 
motion integrated once ( ()5-l() and H5-2() ) are kept unchanged since they are merely 
algebraic relations between and P. Hence, if the evolution of the "constants of 
motion" P is given, we do not need further information about the self-force in order to 
evolve (r(A), 9{X),t{\), (fiX))- Here we give a prescription how to integrate the orbit 
for a long period of time. The properties of orbits in Kerr spacetime which we have 
already mentioned largely simplify the problem of solving the orbital evolution. We 
will find that the leading order approximation can be obtained by using the adiabatic 
approximation as expected. We will also find that the leading order corrections to 
the adiabatic approximation coming from the instantaneous self-force in the linear 
order are given by a few time-averaged quantities constructed from the self-force. 
The basic idea of the discussion in this subsection is found in Ref. 33). 

We begin with discussing orbits for fixed constants of motion P. It is convenient 
to introduce phase functions x*" and by x"' = ^a{P)X, where f2r{P) and Oq^P) are 
angular frequencies of oscillations in r and 9 directions respectively. Note that we can 
choose the initial values of the phases x"(0) arbitrarily. We define functions f{P,x^) 
and 9{P, x^) as solutions of the r-and ^-components of the geodesic equations. Here 
we fix the ambiguity in the choice of phases of these functions so as to satisfy x*^ = 
mod 27r for r{x^) = r_ and = mod 27r for 9{x^) = 9^. Here r_ and 9- are 
minima of r and 9 for given /*. Since the evolutions of r and 9 for < x^ < 
and for vr < x" < ^vr are symmetric, we automatically have x^ = mod 27r for 
r{x^) = r+ and x^ = mod 27r for 9{x^) = 9^ with this choice of phases. Further 
we redefine the functions t^^^ and (p^"^^ here as functions of /* and x^; and similarly 
t(^) and (f'^^^ as functions of P and x^ ■ Again, for definiteness we fix these functions 
so that t(")(r,x") = and (p^'^'^ {P , x") = for x'' = mod 2tt. 

Now we are ready to introduce our parametrization to describe orbits when we 
take into account the self-force. Our proposal is to promote /* and x'^ to functions 
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of A as 

r(A)=f(r(A),x'-(A)), 

e{x) = e{r{x),x'{x)), 

t(A)=tM(r(A),x^(A))+tW(r(A),x'(A))+x*(A), 
V9(A)=(/,M(r(A),x'"(A)) + v^(^)(r(A),x'(A)) + F(A). (5-17) 

Here we have also introduced x^{X) instead of the integrals such as / (dt/dX) {P{X))dX. 
Below we derive equations for P{X), x"(A) and x^{X). 

First we examine the evolution equation for r. As we have mentioned earlier, 
Eqs. (|5-H) and (|5-2|) ) are kept unchanged even if the self-force is taken into account. 
Therefore we have 

±.^?timxn„^^rwi ,5.18) 

since f is the solution of the geodesic equation for fixed On the other hand, taking 
the A-derivative of r(A) in the form given in Eq. (jS-lTf) . we obtain 

Comparing these two expressions, we obtain an equation for dx^ /dX as 



where 



dx''\ ( df \ ^ df dl 



dX J ydx' J dP dX' ^^'^^^ 

Near the turning points r = r±, f can be expanded as r = r±[P) + 0{{Ax^)'^), where 
Ax^ is the difference of x from its value at r = r-i-. Then 6 {dx'" /dX) looks singular 
since dr/dx^ behaves like ~ Ax near the turning points, but in fact this term is 
not singular. As we shall see immediately below, the other factor {df / dP){dP / dX) 
simultaneously goes to at the turning points as long as the self-force stays finite. 

We shall show {df/dP){dP/dX) = 0{{Ax^)'^)- Differentiating the equation 
{dr/dX)^ = R{P,r) with respect to A, we have 

'dr\ d'^r _ dRdP dRdr 

dx) dX^ ~ dP~dX^~d^dX' ^^'^^^ 

Hence, we conclude that {dR{P ,r) /dP){dP /dX) = for r = r±. On the other hand, 
by the definition of r±, we have R{P ,r±{P)) = 0. Differentiating this identity with 
respect to P, and contracting it with dP/dX, we have 

{dR\ Or^dP (dR\ dP 

Hence, except for circular orbits, in which (dR/dr)^^^^ = 0, we establish {dr±/dP) 
{dP /dX) = 0. When the orbit is circular, r is given as a function of P. Hence, we 
do not have to care about the evolution of x^ ■ 
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Now we return to Eq. (|5-2n|) . To evaluate dP /dX, here we accept the use of the 
geodesic momentarily tangential to the orbit in evaluating the self-force. The errors 
caused by this approximation need future investigation. Under this approximation, 
the change of the "constants of motion" becomes a double periodic function in A, 
and is expanded as 

^ = (^) + E ^:.,n«U^(A))e^(--^(^)-^--^(^« + O(M^), (5-24) 

where Jnr,ne{P W) are coefficients to be computed from the instantaneous self-force. 
The last term of 0{fj?) is the correction due to the second order self- force, which 
currently we do not know how to evaluate. In what follows we will give a rough 
estimate of the error due to this second order contribution, but its derivation is not 
rigorous at all. The second term in (|5-2n|) is also a regular double periodic function. 
Hence it accept a similar expansion. We can write Eq. (|5-2()j) as 

+ Yl '5I?"'-'"''(/-''(A))e^("'-^''(^)+"«'^'(^)) +0(^2)^ (5.25) 

{n,.,ne)7^{0,0) 

The effect of the third term looks as large as the second one at first sight. However, 
if we integrate the above equation once, one finds that the second time grows as 
0{nAX), where AX is the length of integration time. Since we are interested in long 
term orbital evolution, AX is assumed to be large. On the other hand, the third 
term can be integrated by parts as 



i{nr{dx''/dX) + neidx^dX)) 

J dX dlH{nr{dx'-/dX) + ne){dxVdXy ^ ' 

The first term in the last expression does not grow even for a large value of Z\A, 
while the second term becomes higher order in ^ due to the appearance of the factor 
dP/dX. We express the result of integration of (|5-25j) as 

X'-(A) = I |f?, + ^5(^^)^|ciA + O(M^A)0) + O(M'Z\A). (5-27) 

The last term is the correction coming from the second order self-force. Discussion 
about the ^-component is completely parallel to the r-component. 

Next we consider the t-component. By the definitions of t^^'^ and t^^\ we have 

dt atM^ /dt\ , 



dX dy dx^ " \dX 
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On the other hand, differentiating t{X) in the form of (|5-17|) . we obtain 
dt _ dP dx'' dP dt^^^ dx^ dx* 

From a comparison of these two equations, we find that the evolution of is to be 
determined by 



(7 ~t\ ^r,ng 



(nr,ne)^(0,0) 



with 



y 91* rfA Sx'- [dXj dP dX dx' [dXj- ^^''^ ' 

Then, this equation can be integrated as in the case of x"- A completely parallel 
discussion goes through for the (/j-component, too. 

Now we consider the evolution of P{X). The evolution equations for /*(A) are 
already given by Eq. H5-24() . We express the result of integration of (|5-24|) as 

P{X) = P^,{X)+SP{X). (5-32) 

The first term Pj^^{X) is the contribution only from the first term on the right hand 
side of Eq. (|5-24|) . We find that the contribution to 5P{X) from the terms including 
Jnr,ne remains 0{n{AX)^) or higher. On the other hand, the second order self-force 
contribute to 5P as terms of Odi^AX). 

The errors in /*(A) propagate to x"(A) and x'^(A). By using the same argument 
that we have already used many times, the propagated errors caused by the terms 
which include Jn^^ng most 0{ij,{AX)'^), while the effect of second order self-force 

can be as large as 0(/i^Z\A^) for both x°'W and x^{X). To conclude, we found 

r (A) = /^,(A) + 0{f,{AXf) + 0(/.MA), 

x"(A) = xUA) + / (^S (^)^dA + O(M^A)0) + O(/iMA2),j 



with 



X^(A) = xl(A) + 1(^5 yx + O(^(Z\A)0) + Oif^'AX'), (5-33) 

Ia,W ■■= J (^)(^Ad(A))dA, 
X^,(A) := I Qa{pAAm\ 

^'^'^^^ I {^) ^Ad(A) := / (^) (Ad(A))«!A. (5-34) 
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The last term for each /*(A), x°'{^) comes from the second order self- 

force. The second terms in the expressions for x°'W a.nd x^{X) are 0{fiAX). When 
fiAX <C 1, the deviation from the adiabatic approximation is small. In this sense 
adiabatic approximation is a good approximation. On the other hand, once fiAX 
becomes 0(1), we cannot neglect the second order self-force, either. Hence, roughly 
speaking, the orbital evolution which takes into account the instantaneous self-force 
but only at the linear order will not be a better approximation than that obtained 
by using the adiabatic approximation. Even when fj,AX is moderately small, the 
leading corrections are given by quantities averaged over a long period of time, 
{S (dx^/dX)) and (5 [dx^ / dX)) . Again we can repeat the same argument that was 
used to justify replacing the retarded field with the radiative one when we evaluate 
dQ/dt. Therefore adiabatic approximation is sufficient to calculate these averaged 
quantities. 

Nevertheless, we are not saying that the study of the instantaneous self- force 
is not important. When we consider the whole process of inspiral of binaries, the 
time scale is inversely proportional to /i. Therefore, both the terms of 0(/iZ\A) 
and of 0{fj?AX'^) potentially cause measurable effects. To evaluate these effects, we 
need to know (maybe some time average of) the second order self-force, and for this 
purpose the study of the instantaneous self-force at linear order will be necessary. 
Furthermore, there might be situations in which the effects of the oscillating part 
of the self- force, i.e., the terms which contain Jn^^ng ^-^^^ significantly enhanced. In 
the above discussion, when we perform integration by parts such as shown explicitly 
in (|5-26j) . a factor nr{dx^ /dX) + n0{dx^ / dX) appeared in the denominator. As the 
"constants of motion" P evolve, this denominator eventually can cross for a certain 
combination of and uq. In such a case, the contribution from Jnr,ng ™ay leave some 
significant effect on the orbital evolution. '^^^ Suppose that nr{dx^ / dX) + ng{dx^ / dX) 
vanishes at A = Aq. Such a stationary phase point (SPP) will cause an additional 
shift in /*(A) which will be estimated by the Gaussian integral 

Alir,T. ^ I dXf fji(';^))eiKx'-(Ao)+nex'(Ao))+f{n,(d2xV'^A2)+ne{d2x''/rfA2)}(A-Ao)2 



~W TTT^TT^ ,n ,,t J'n n (P(Ao))e^('^'-'^''(^o)+"''>^*(^«)). (5-35) 

V dp d\ '''S dp d\ ) 

This correction is 0{y/Jl), and it will induce shifts in phases and x^ of 0{^/JlAX). 
Hence, this naive order counting indicates that this correction is the leading order 
correction to the adiabatic approxiamtion. However, when we consider relatively 
non-relativistic orbits with small eccentricity, such stationary points will appear only 
when the values of Ur and/or ng are large. In such cases the coefficients Jnr,ng ^i^^ 
be suppressed by some large powers of the eccetricity or (v/c). Hence, the effect 
will practically remain small. To the contrary, when we consider highly eccentric or 
highly relativistic orbits, the corrections due to stationary points may become really 
0(-^/7tZ\A) without any significant additional suppression. Again, to evaluate these 
effect quantitatively, we need to compute the instantaneous self-force. 

Before closing this subsection, we want to emphasize that the errors in phases 
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and can grow in a dynamical time scale if we do not use an appropriate integration 
scheme as presented here. Namely, they can be as large as 0(/x(e"'^'^ — 1)) with a 
being a constant of 0(1). 

§6. Toward post-Teukolsky formalism 

Now we consider to extend the Teukolsky formalism to the second order, which 
seems to be indispensable to improve predictions for waveforms beyond the level of 
adiabatic approximation. Perturbed Einstein equations up to the second order take 
the form 



h(^^ =Tl^-G^ h«,h« , (6-1 



where is nothing but GjfJ, in Eq. (|3-5jl . Projection of this equation can be done 
formally as in the case of linear perturbation as 

LC(2) = ^r(2). (6-2) 

Here ^^"^^ is defined in the same way as C^^) just substituting h^^) with h(^). The 
second order source term T*^^-* has a spatially extended distribution due to the non- 
linearity of gravity. The differential operator L is the same as the one for the linear 
perturbation. Therefore all the technical difficulties which arise in the second order 
for the first time are in how to evaluate the source term. To obtain the corrections to 

(2) 

the energy momentum tensor T^u , we need to know the correction to the trajectory 
of the particle taking into account the self-force. Hence as a first step toward the 
post-Teukolsky formalism, we examine the self-force. 

6.1. gravitational self-force 

When we consider the point particle limit, the full self- force diverges at the 
location of the particle, and hence needs to be regularized. It is known that the 
properly regularized self- force is given by the tail part of the self-field, which is 
obtained by subtracting the direct part from the full field. Here we do not give 
the precise definition of the direct part, but, roughly speaking, it is the part that 
propagates along the light cone unaffected by curvature scattering. Since direct part 
does not contain curvature scattering effect, the direct part of the local field near 
the particle is solely determined by the local geometry and the orbital elements of 
the particle. The justification of this prescription is given in^^^ for the scalar and 
electro-magnetic cases, and in^^)'^^) for the gravitational case. An equivalent but 
more elegant decomposition of the Green function was proposed by Detweiler and 
Whiting, ^^^'^^^ in which the direct part is replaced by the S'-part and the tail part 
by the i?-part. The S'-part is defined so as to vanish when two arguments x and x' 
are timelike. When S'-part is subtracted from the full field, the remainder i?-part 
gives the regularized self-force. The advantage of this new decomposition is that the 
S-part is symmetric with respect to x and x', and it satisfies the same equation as 
the retarded Green function does. This implies that the i2-part now satisfies the 
source-free, homogeneous equation. 
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6.2. subtraction and regularization^^^ 

Since we do not know a direct way to compute the i?-part, we instead compute 
[ip^] (r) = F" [V"^""] (r) - [^^] (r). Since both the terms on the right hand 
side are divergent, this expression does not make sense unless we regularize the di- 
vergent quantities. A practical way of regularization is the so-called "point splitting" 
regularization. We evaluate the expression for the force not exactly at the location 
of the particle but slightly off the point. If we subtract the /S-part appropriately, the 
coincidence limit must be well-defined. 

For an actual computation, we decompose divergent full- and S-parts of the force 
into terms labelled by the total angular momentum ^^37)-45) 

oo 

F« (r) = hm V (Ft U^""] (r) - F," [V^^] (r)) . (6-3) 

Then each ^-order term stays finite even in the coincidence limit. Here the question 
is whether we can change the order of the two operations as 



^x^.(r)V L J J 

In general these two operations do not commute. 

For example, let us define two functions A{x) and B{x) as A{x) = Yle^oi^ ~ 
and B{x) = Yle^ii^ ~ fo^" x>i). For a; 7^ 0, both A{x) and B{x) are convergent. 
Therefore we can compute the expression A{x) — B{x) unambiguously, and it is 1. 
In this case it is trivial that the difference A{x) — B{x) = 1 -|- X^^x *-* ^® uniformly 
convergent, and each £-th order term is continuous. Therefore, one can change 
the order of two operations. Irrespectively of the order of operations we arrive at 
the same answer. The importance of the conditions of uniform convergence will 
become clearer if we consider the case with B{x) = X^^i(l — xY~^. In this case 
A{x) — B{x) = 0. Nevertheless, if we take the limit x — > first, we end up with 
A(0)-S(0) = 1. Now the term at the Mh order is {l-xf -{l-xf-'^ = -x(l-a;)^-^ 
As X becomes closer to 0, the convergence becomes worse. This series is convergent 
at each point x but not uniformly. 

A practical way to guarantee the uniform convergence in the harmonic expansion 
of the force is to use the same prescription for both /w/Z-part and S'-part to extend 
the force off the trajectory and to use the same harmonics. The S-part of the force 
is determined by the local expansion near the particle. It is composed of terms like 

Y.^^f-^cd[z''{r\u>^{T)l (6-5) 

Here, the separation from the trajectory of the particle {x — z{t))^ was denoted by 
{T,R,O,0), and e is the spatial distance between x^ and the trajectory. This type 
of function can be expanded in terms of spherical harmonics.^^^'^^^'^^) On the other 
hand, the full-part composed of the master variable is also computed by using the 
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spherical harmonic decomposition*-* . After averaging over the angular direction from 
the location of the particle, we will obtain the full- and the S-parts of the force in 
the form of 

oo 

F-{t,R,X) = J2Ft{r,R)Pe{cosX), (6-6) 

e=o 

where X is the angle between and z^{t). These expansions have a little ambiguity 
since the extension of the self-force is given only locally, in the vicinity of the particle. 
Nevertheless, the ambiguity does not affect the force evaluated near the particle after 
summation over £ since that is the basic requirement for the construction of harmonic 

expansion."^ ''^ ' "^^^ ' '^"'^^ 

Now we consider the difference between the full-pait and the S'-part of the self- 
force. If both parts are computed appropriately, the difference must become i?-part, 
which is regular at least near the location of the particle. Therefore, the series 

oo 

E i^fuuAr, R) - Fi^T, R)) P,(cos X), (6-7) 

e=o 

must be uniformly convergent. Thus we can take the coincidence limit, and we arrive 
at 

oo 

nir) = Yl i^fullA^^ = 0) - F^^,{t, R = 0)). (6-8) 

e=o 

The S'-part and hence also full-part of the self-force in the harmonic gauge takes 
the form of37)>39),4i),45) 

1\ „„, 1^-^ 



F^{T,R = 0) = A"\^£+-j+B" + C"\^£ + -j -\-Df, (6-9) 

The constants A", i?", C" and the residual part := J2e^o called "regular- 

ization parameters" . The choice of the residual for each £ mode, Df, can be changed 
rather freely by changing the behavior at a large separation angle, but the value of 
regularization parameter I?" is not affected. The first three regularization parame- 
ters A", B" and must be common for both full- and S'-parts since otherwise the 
i?-part diverges. Hence, what we need to evaluate is the difference 

mr) = Y.i^%ii,e-^s,e)- (6-10) 

e 

The terminology "mode-sum regularization" or "mode decomposition regular- 
ization" is often used for the above prescription. However, we want to stress that 
the mode decomposition itself is just a useful tool to simplify the procedure of the 
point splitting regularization. 



*^ The way how to extend the force must be carefully chosen so that the harmonic decomposition 
of the full-part of the force is computable. Alternative way to avoid this problem is to subtract 
S-part at the level of the metric perturbations®^^ (or at the level of the master variable''®^). 
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6.3. gauge problem and use of intermediate gauge. 

The formal expression for the self- force is derived in harmonic gauge, and S- 
part is also given by using a local Hadamard expansion of the Green function in the 
harmonic gauge. However, when we calculate the full metric perturbations by using 
some metric reconstruction method from a master variable, the gauge that we can 
practically use is restricted to Regge- Wheeler- Zerilli gauge for Schwarzschild case or 
radiation gauge for Kerr case. Since the force is a gauge dependent quantity, regu- 
larization parameters can be different in different gauges. Hence, simple subtraction 
of 5-part for harmonic gauge does not work. 

Here we present the basic idea of the intermediate gauge method^''^ in a manner 
slightly modified from the original. We associate subscripts, full, R, S to denote full 
perturbations, R-pavt and S'-part, respectively. In any gauge G we can define the 
i?-part as 

:=hg/-hf . (6-11) 
We consider a gauge transformation from the harmonic gauge to the gauge G: 

h(^) = h(^) + VC(^-^^)[h(^)], (6-12) 

where (H) represents harmonic gauge, and denotes the change of metric gener- 
ated by an infinitesimal coordinate transformation ^ — ^(^^'^)[h(^)] is 
the generator that transforms the metric perturbation in harmonic gauge h^^^ into 
that in the gauge G. Then i?-part of metric perturbations in the harmonic gauge 
can be rewritten as 

hSf^=hg,-hf)-v4^-^\ (6-13) 



where 
and 



hf)=h(f)+v4^^^\ (6-14) 

Here ^^^^'^^ is not specified yet, and it is not necessarily equal to ^'^■'^"^'^^[h^'']. 

(C) (G) 

This is the reason why we used the notation instead of . 

The idea is to drop the last gauge term V^^^^^^ [h^''] in H6-13() . Namely, the 
force is computed from the regularized metric 

tr=42,-hr'. (6.16) 



A trajectory in this intermediate gauge is related to that in the harmonic gauge via 
the gauge transformation specified by C/?^"^^^- For the intermediate gauge to be 
useful, it should satisfy the following two conditions. 

1. We need to choose ^5^^'^^ so that h^^ ^ completely cancels the singular part in 
h^^]^ in the coincidence limit. 

2. There is no secular growth in the generator of the gauge transformation 
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The condition (1) is the requirement that is a good approximation of the 

true S'-part gauge transformation, ^'^■^^'^^ [h^'']. If the difference between i^^^^^ 

and ^(■'^~*'^)[h^^] is regular and finite, we have = [h^^] except for 

regular term. Hence, the coincidence limit of is guaranteed to be finite. To 

obtain a sufficiently good approximation of ^(^~**^^[h^^''], we just need to know h^^"* 
as a local expansion near the trajectory, and that is the best we can do. 

Here it is very convenient if the gauge G satisfies a property that the gauge pa- 
rameters to transform metric perturbations in another gauge to the specified gauge 
G are determined without temporal or radial integration. This means that the equa- 
tions to determine the gauge parameters are solved locally on a sphere. Then trans- 
formation into G-gauge is fixed unambiguously just by looking at the perturbations 
in the vicinity of a sphere. We call such a gauge as a Gauge Operationally Deter- 
ministic on a Sphere (GODS). Regge- Wheeler gauge is GODS.^^^ Radiation gauge 
in Kerr case does not seem to be GODS in its original form, but there seems to be a 
modification which allows it to transform into GODS in the sense of expansion with 
respect to Kerr parameter a*\ 

In understanding the true difficulty about the gravitational radiation reaction, 
it is very important to realize that the condition (2) is really necessary. First of all, 
if the condition (2) is not required, we can add any finite gauge transformation. As 
a result, the trajectory in the coordinate representation can be arbitrarily changed, 
which is quite unsatisfactory situation. If the amplitude of gauge parameters is 
large, the coordinate values of the particle in H and G gauges are quite different. 
Even in that case, although it is quite counter-intuitive, one can cancel the divergent 
pieces in metric perturbations as far as linear perturbation is concerned, in which the 
gauge transformation does not couple with metric perturbations. However, once we 
consider the second order perturbation, large gauge parameters will cause disaster 
to the perturbative expansion. 

Now we come back to the issue how to guarantee the condition (2). When 
the gauge G is GODS, we can define ^5^"^*^^ in terms of local quantities without 
including any integration over i or r. With such a choice of , we do not 

have to worry about the condition (2) for 5-part since there is no secular growth in 
Cs What we need to guarantee to satisfy the condition (2) is the absence of 

secular growth in ^^^^'^^[^^f^ii]- 

We explain one practical scheme to guarantee the absence of secular growth in 
^(if-+G) [h^^^-'j**). We first recall that hj^^-*^ must stay finite by assumption. If this as- 
sumption does not hold, the whole formulation of self-force based on harmonic gauge 



*^ We would like to come back to this point in future publication. 

There is another way given by Amos Ori during post-Capra discussion meeting (2003 Kyoto), 
which does not rely on the presence of GODS. We give a brief explanation of his argument here 
although it might be inaccurate. When the source can bo decomposed into Fourier mode in time 
direction, both h^^j'; and h^^], corresponding to a partial wave will be periodic. Namely, we assume 
h/j,(((t,x) = hfuii{t + T, x). Then, the gauge transformation connecting between two such metrics 
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breaks down, and hence we lose the whole foundation. For GODS, ^'•^^'^^ [hj^^J is 
guaranteed to stay finite since so ii^f^ji is. 

§7. Other topics 

Here we briefly mention a few topics which we have not yet discussed at all. 

7.1. analytic approach to the self-force 

There are a few calculation of the instantaneous self-force using numerical ap- 
proach for a limited class of orbits. If we do not use the fully numerical approach, 
the full-part of the self-force is given in the form of Fourier expansion in time. On the 
other hand, the 5-part is given in the local expansion, and hence is expressed in the 
time domain. To perform the subtraction, we need to transform the full-part (or the 
S-part) into the expression into the time domain (or into the frequency domain). As 
was mentioned earlier, a systematic method to solve the linear perturbation equation 
analytically is already known. ^^^'^^^ Taking advantage of it, we proposed a method 
to perform this transformation between different domains analytically. 

7.2. £ = and £ = 1 modes 

There was a debate on how to treat £ = and £ = 1 modes. ^''^'^^^ Here we give 
a simple-minded understanding of this complicated issue from the view-point of the 
intermediate gauge approach, although deeper understanding might indeed become 
indispensable when we consider the second order perturbation. 

The treatment of these lower lying modes in black hole perturbation theory is 
different from the other modes. The Regge-Wheeler-Zerilli formalism can handle 
these modes separately, but their equations are not hyperbolic. Therefore the re- 
tarded boundary conditions do not fix the boundary conditions for these modes. In 
the Teukolsky formalism, these modes are absent from the beginning. 

Those £ = and 1 modes are composed of physical part and gauge part. In the 
intermediate gauge approach, no debate can arise in the gauge part. As long as it 
does not secularly increase, we do not care about it at all. The problem arises only in 
the point how to determine the physical part. There exists some information which 
is not encoded in the master variable, which is what we call £ = and 1 modes. 

The particle motion defines a three dimensional tube in four dimensional space- 
time specified by the trajectory {t,r) = (tz{T),rz{T)). This tube divide the back- 
ground spacetime into two pieces: one containing the infinity and the other con- 
taining black hole horizon. 

In both inner and outer regions, the metric perturbation is a homogeneous so- 
lution of Einstein equations. It is given by the part reconstructed from the master 

V(^^i,) satisfies the same periodicity. Hence, we have [(,^-.1^ + x) = [^^i-i, + ^i^:tj.]{t + T, x). This 

implies that 5^(t-|-T, x) — {^(t, x) satisfies the Killing equation. Therefore, — (periodic piece) + 
tK^, where is a Killing vector. For black hole background the variatiou of is limited, and 
hence for any choice of non- vanishing K^, the gauge transformation V(p(t-ft'i,)) docs not vanish at 
r* — > ±oo. As far as the h^^], is guaranteed to go to zero at infinity or on the horizon, we find 
= 0. Hence, the /itZZ-part of the gauge transformation does not have any secular growth. 
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variables with i >2 with an additional piece which does not affect the master vari- 
ables with £ > 2. Possibly such additional perturbations are those which are given 
by a change of the parameters contained in the background metric: 



Q BH Q BH 

''^'^ " dM + da 



Sh^^ = -T^SM + ^7f-5a, (7-1) 



where M and a are the mass and Kerr parameter of the central black hole, respec- 
tively. For vanishing of the master variable, no other perturbations are possible. The 
Schwarzschild case should be considered as a special case of Kerr. If we consider an 
orbit with inclination, one may think that these two parameters are insufficient be- 
cause we do not have parameters corresponding to the rotation in x and y axes. 
However, these rotations other than that in the z-direction can be absorbed by a 
gauge transformation, i.e., by global rotation of angular coordinates. 

If we do not care about small error of 0(/x) in the estimate of the mass and the 
angular momentum of the central black hole, SM and 6a are unimportant. They are 
renormalizcd in the definition of the mass and the angular momentum of the central 
black hole. Hence, this issue about i = and 1 is not an issue of debate at the 
lowest order in ^u. It becomes an issue only when we discuss it in connection with 
the second order perturbation. 



§8. Summary 



In this paper, we gave a brief review of the recent development in the study 
of gravitational radiation reaction problem in the context of generating gravitation 
wave templates, with some new insights. 

First, we reported the adiabatic approximation to the radiation reaction, in 
which long time average is assumed in evaluating the change rates of the "constants 
of motion" . The radiation reaction to the Carter constant had been a long-standing 
issue, but now we are ready to compute it in the adiabatic approximation. We 
explained the formulation to evaluate the change rate of the Carter constant. We 
have presented a method to integrate the evolution of orbits for a long period taking 
into account the radiation reaction. We have shown that the errors in the phases 
of the orbits obtained by using the adiabatic approximation are 0{iiA\, //^Z^A^) or 
higher, where fi is the mass of the small compact star orbiting the central black 
hole and AX is the time duration of integrating the orbit. In this order counting, 
AX is supposed to be large, but iJ,AX is small. Typically, the AX oc //"^ since the 
evolution due to radiation reaction is slower for a smaller mass. The phase errors 
caused by ignorance of the instantaneous leading order self-force is 0{iiAX). But 
when those errors become large, the contribution from the second order self-force of 
0(/Li^Z\A^) becomes comparable. However, there is a possibility that further study 
on the leading order self-force alone, beyond the level of adiabatic approximation, 
can improve the gravitational wave templates dramatically in some cases. As the 
orbit evolves adiabatically, the frequencies of oscillating part of the self-force change, 
and eventually one of them may cross zero. A rough order of magnitude estimate 
suggests that the corrections in phases due to this accidental appearance of zero 



Gravitational radiation reaction 



25 



frequency modes can be as large as 0(y7^Z\A). For more definite estimate of this 
effect we need to study the instantaneous self-force. 

Furthermore, once we start to detect gravitational waves from binary inspiral 
with an extreme mass ratio, the detailed comparison of the signal with the theoretical 
prediction will become possible. In such a situation, we will wish to have a theoretical 
tool which can predict the phase evolution with an accuracy less than 0(1). For this 
purpose, we need to understand the self-force not only at the leading order but also 
up to the second order, although probably some kind of averaged values will be 
sufficient for the second order self- force. 

Also as a step toward this goal, complete understanding about the first order 
self-force will be necessary. Recently, there have been a lot of developments also in 
this direction. We have presented here a biased summary about this issue, partly due 
to my understanding and also due to lack of space. So far, the lowest order self-force 
has not been discussed extensively as a step toward the second order self-force. ^'^^ In 
order to reflect the stored knowledge about the self-force to the improvement of the 
theoretical prediction of gravitational waveforms, we need to develop a formalism to 
evaluate the second order self-force. 
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